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I. INTRODUCTION 

The two dimensional (2D) fully frustrated XY model 
(FFXY)i42, is one of the most intriguing of "simple" 
statistical mechanics models. The doubly degenerate 
checkerboard pattern of vortices in the ground state 
leads to an Ising-like discrete Z(2) symmetry in addi- 
tion to the Kosterliz-Thouless-like continuous 0(2) sym- 
metry associated with the uniform rotation of all phase 
angles^. It remains controversial whether there are two 
distinct phase transitions Tkt < Tj, with Tkt mark- 
ing the breaking of the 0(2) symmetry and Tj mark- 
ing the breaking of the Z(2) symmetrji^i^, or rather a 
single transition in which both symmetries are broken 
simultaneously . 

Recently, Korshunov^ presented an argument for a new 
interfacial transition T^^ in the 2D FFXY model, lying 
well below the bulk transition(s), arising from the un- 
binding of step excitations of unit height ( "kink-antikink 
pairs" ) on the domain walls associated with the Z(2) sym- 
metry. Korshunov argued that the kink-antikink unbind- 
ing transition leads to a decoupling of phase coherence 
across domain boundaries, supporting the identification 
of the FFXY with the coupled XY-Ising modelS. Kor- 
shunov further argued that this effect necessarily leads to 
the scenario of two distinct bulk transitions, Tkt < Tj. 

Earlier, Lee and co-workers, first in simulations of the 
2D FFXY model with Langevin dynamicsiS, and then 
in its dual Coulomb gas (CG) model with Monte Carlo 
(MC) dynamicsii, found evidence for a transition in do- 
main wall morphology in simulations of the ordering ki- 
netics of domain growth following a sudden quench. They 
interpreted this as a finite temperature roughening tran- 
sition of the Ising-like domain walls. Jeon et al}^ made 
similar conclusions in simulations of the 2D FFXY with 
resistively shunted junction dynamics. Korshuno\i^ how- 
ever, has argued that domain walls should be rough at 
all temperatures. 

In this paper we present the first direct numerical ev- 
idence demonstrating the existence of the kink-antikink 
unbinding transition at a temperature T-^j below the bulk 
transition(s). In agreement with Korshunov's predic- 
tions, we show that phase angles on opposite sides of the 



domain wall decouple above T^. The numerical value we 
find for Tu, is comparable to that of the morphological 
transition found by Lee et al^, however we explicitly 
demonstrate that domain walls are rough at tempera- 
tures well below T^,. This indicates that the transition 
seen by Lee and co-workers is really kink-antikink un- 
binding, rather than roughening. 



II. THE MODEL 
A. The fully frustrated XY model 

We study the 2D FFXY model on a square lattice, 
given by the Hamiltonian^'^, 



H[e(vi)] - ^ y (0(r, -I- A) - 6{v,) - A,{y{)) 



(1) 



Here 9{Yi) is the thermally fluctuating phase angle of the 
planar XY spin on site = XiX -\- yiy {xi, yi integers) 
of a La; X Ly periodic square lattice, pL = x,y label the 
bond directions of the lattice, and ^^(r^) is the quenched 
gauge field on the bond leaving site in direction /t 
(with A-fj^{Yi + jl) = —Af^{ri)). For full frustration, the 
A^{ri) are constrained so that their directed sum going 
counterclockwise around any plaquette P of the lattice is 
fixed (modulus 2tt) to. 



^^^(r,) = ^ 



(2) 



To implement the constraint of Eq. Q , we use the spe- 
cific gauge choice, 



A,(r,) = 0, A,(rO = (-!)"' (7r/2) 



(3) 



The interaction potential V{4>) is periodic on [0,27r), 
with a single quadratic minimum at = 0. We will take 
for V{(j>) the commonly used Villain functioni^. 



V{(f>) = -Tin 



(4) 
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FIG. 1: Various configurations of tlie domain wall in a L x 
{L + 1) system: (a) ground state, (b) finite widtli step of unit 
lieiglit (kink-antilcink pair), (c) isolated kink of unit height, 
(d) isolated kink of height two. A (+) indicates the presence 
of a vortex in the XY model, or a charge qi — 1/2 in the dual 
Coulomb gas; a ( — ) indicates the absence of a vortex in the 
XY model, or a charge g; = —1/2 in the dual CG. x is the 
horizontal direction, and y is the vertical direction. 




FIG. 2: Variation of total free energy F with total twist Ay 
applied transverse to the Ising-like domain wall, for two dif- 
ferent temperatures in a system of size L = 128. 



The boundary conditions for the phase angles are, in the 
most general case, given by^iii^ 



(rj + Lf^fi) - 0(ri 



(5) 



where £ [0, 27r) is the total twist applied across the 
system in direction jl. — corresponds to peri- 

odic boundary conditions. Alternatively, if one makes 
the change of variables, 0'{ri) = 9{ri) — r.^ • d, with 
= A^j/Lp, then the system has periodic boundary 
conditions in the 9'{ri) and the applied twist appears 
as an additive constant to the gauge field, ^^(ri) 

To study the behavior of the Ising-like domain walls 
we consider systems with sizes = L, Ly — L -\- I, 
with L even. The odd length Ly forces into the ground 
state checkerboard pattern of vortices a single straight 
domain wall running the length of the system in the x 
direction. This is illustrated in Fig. where a {+) 
signifies a vortex in the phase angles 6{ri), and a (— ) 
signifies the absence of a vortex. 

Phase coherence in the FFXY model is most conve- 
niently studied by considering the dependence of the total 
free energy F on the total twist A^ applied across the sys- 
tem (see Eq. Q). In a phase coherent ordered state, we 
expect that F{A^) varies with the twist A^; in a phase 
incoherent disordered state, we expect -F'(A^) is indepen- 
dent of A^ in the thermodynamic limit of L oo. The 
dependence of the free energy on A^ is readily obtained 
by using fluctuating twist boundary conditioned'^ , in which 
one treats the applied twist A^ as a thermally fluctuat- 
ing degree of freedom. If Z is the partition function for 
this ensemble, then the probability P(A^) of finding a 
state with a particular twist A^j is given by. 



(6) 



and so the free energy with respect to a reference twist 
is, 

F(A^)-F(A,,o) = -Tln[P(A^)/P(A,,„)] . (7) 



A^o is, 



The probability P(A^t) is directly measured within our 
fluctuating twist Monte Carlo simulation. We choose the 
reference twist A^o to be the value of the twist that min- 
imizes the free energy P(A^). For the gauge choice of 
Eq. |j3J), it is straightforward to see to that the minimiz- 
ing twist in the x direction is at Axo = 0. In our simula- 
tions we keep a fixed twist A^; = 0, and consider only the 
dependence of the free energy on the varying twist Aj,, 
transverse to the Ising-like domain wall that is introduced 
in our L x (L+l) systems (see Fig.^). In Fig.|21we show 
sample results from our simulations for F{Ay) — F{0) vs. 
Ay at two different values of P < T^, for a system of size 
L — 128. We see that P(Aj,) has two equal minima at 
Aj^ — and tt (i.e. periodic and antiperiodic boundary 
conditions). One of these minima corresponds to states 
where the domain wall sits at even values of the height y, 
while the other corresponds to states where the domain 
wall sits at odd values of the height y. Below the kink- 
antikink unbinding transition , the system is in a state 
of broken translational symmetry; due to the free energy 
barrier between the two minima, states in which the do- 
main wall is at an even height cannot be reached from 
states in which the domain wall is at an odd height. As 
noted by Korshunov^, this broken symmetry is restored 
when phase coherence transverse to the wall is lost, i.e. 
when F{Ay) becomes independent of Aj,, and so the free 
energy barrier between Ay = and Ay = tt vanishes. 
Alternatively viewed, when the domain wall changes its 
height by an odd number, the system acquires an average 
twist of TT in the y direction. Thus, restoring the sym- 
metry of domain wall translations leads to phase angle 
fluctuations that destroy phase coherence transverse to 
the direction of the wall. 

In our numerical work we will use two convenient mea- 
sures of the variation of P(A^) with A^. The first is the 
helicity modulus^ii, T^, which measures the curvature 
of P(Ap) at its minimum. 



T^(Px, Ly 
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where = 4>^iri) = 0{ri + fi) - 6{ri) - ^^(ri), V 
and V" are the first and second derivatives of the Villain 
function of Eq. and (. . .)o indicates a thermodynamic 
average in the ensemble with fixed twist = 0. A 
second measure is, 



F{it/2)-F{Q) 



values of i^(Ay), as Aj, is varied at fixed A^; = 0. Since 
is an intensive quantity, it should approach a value 
independent of system size as L — > oo. The parameter 
/S.F scales as T{A/ L)'^L^ in D dimensions, and so for 
_D = 2 it also becomes independent of system size as 
L oo. 



B. The Coulomb gas 

Although our simulations are carried out in the XY 
variables 0{ri), it is helpful to consider the situation from 
the viewpoint of the dual CG model of logarithmically 
interacting half integer chargesii^. For the case of a fixed 
total twist A^, the XY Hamiltonian of Eq. |^ maps onto, 

ncG = no + Hi . (10) 

Ho is the logarithmic interaction of the charges. 



Ho = ^(27rJ) V qiG{ri - Tj jqj 



where qi — ±1/2 are the half integer charges, neutrality 
is imposed, J^i 1i — 0, and G(r) is the 2D periodic lat- 
tice Coulomb potential^'^^ with G{r) ^ — ln|r| for large 
1 "C |r| ^ L/2. Til arises from the fixed twist boundary 
condition and is given byiii^. 



Hi^VxlA 
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-Vy A, 



Ay 



2'Kpx 



0. However, for this ground state one has Px = i/4, 
hence the ground state energy is minimized when Aj^ = 
[L + l)7r/2 — 7r/2 = Ltt/2. If the location of the domain 
wall was shifted by one unit in height, then the ground 
state would have px = —L/A, and the energy would be 
minimized when Aj, = {L + l)Tr/2 + TT/2 = LTr/2 + n. For 
L even, as we have required, these two values, modulus 
27r, are just equal to and tt. 

The helicity modulus Ty/J maps^ti^ii^ onto the in- 
^9) verse dielectric function of the CG. As noted by 
KorshunoAi^, in order for the domain wall to move a unit 
lattice spacing in height, it must first form a unit step 
of finite width £; £ must be even to preserve charge neu- 
trality (see Fig. ^). We denote the left hand edge of 
the step as the kink, and the right hand edge as the 
antikink. As the kink and antikink separate out to in- 
finity, the domain wall moves one unit in height. As 
shown by Halsey^^, a corner in a domain wall carries 
with it a net charge of ±1/4. The kink, consisting of 
two successive corners with equal +1/4 charge, carries a 
net charge of qkink = +1/2; the antikink carries a net 
charge of gkink = ^1/2- At low temperatures, the loga- 
rithmic attraction between the kink and antikink charges 
keeps them bound with a largest separation ^max(T). At 
higher temperatures, entropy wins out over energy, and 
there is a kink-antikink unbinding transition at , where 
^max(?u,) — > oo. Abovc T^, the kink-antikink unbinding 
leads to diverging dipole fluctuations in the x direction, 
driving (and hence Tj,) to zero. 

The problem of logarithmically interacting charges in 
one dimension (ID) has been treated by Bulgadae\i2&. 
Koshuno\i2i has applied these results to the unbinding 
transition of the kink-antikink pair along the one dimen- 
sional Ising-like domain wall. To include the screening 
effect of charge excitations in the bulk of the system on 
either side of the Ising-like domain wall, we take as the 
coupling between kink-antikink pairs separated at large 
distance to be the helicity modulus of the FFXY for an 
ordinary L x L system, T(L, L), in the limit of large 
enough L. Applying Bulgadaev's exact result for the un- 
binding transition temperature, we conclude. 



(11) 



(12) 

where Vx and Vy are Villain functions as in Eq. but 
with couplings Jx = J(Ly/Lx) and Jy = J{Lx/Ly) re- 
spectively, p is the total dipole moment. 



E 



and Al = T.xM^,y = 0), Al = Y^yM^ 

the gauge choice of Eq. ||2J) we have, 



^1 = 0, A! 



(i + l)7r/2 



(13) 
0, y). For 

(14) 



For the ground state as illustrated in Fig.^, one has 
Py — 0, and so again it is easy to see from Eq. H12|) that 
the total ground state energy is minimized when A^; = 



27rT(L,L)g^,„, 
T 



2, 



-T{L,L) . (15) 



One can reproduce this result using a Kosterlitz- 
Thouless-like argument^ as follows. In analogy with Lee 
et a/. Ai, we consider the total free energy to have a sin- 
gle "free" (i.e. unbound) kink in the domain wall (see 
Fig. ^). Fixing the kink at a given position on the do- 
main wall, its free energy (averaging over all other fluc- 
tuations) is that of an isolated +1/2 vortex in a medium 
with phase stiffness T{L,L); here we use Bulgadaev's 
result^" that kink-antikink pairs in ID do not lead to a 
renormalization of the kink-antikink interaction, and so 
any screening of their interaction is due to charge excita- 
tions in the bulk on either side of the domain wall, and so 
accounted for by the large L value of T. As L ^ oo, the 
leading contribution to this energy is E = 7r(7^jj^j.T InL. 
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FIG. 3: Helicity modulus vs. T for different system sizes 



Lx X Ly. As L 

the same curve. 



oo, Tx{L,L + 1) and Tx{L,L) approach 



The entropy of the kink is just that associated with its 
position along the domain wall, S = —InL. Combining 
gives i^kink = E-TS = {ttT/A-T) InL, which as i ^ cx) 
gives the instabihty temperature for the formation of free 
kinks as = 7i'Tr/4, in agreement with Eq. H15() . 



FIG. 4: Helicity modulus Ty vs. T for different system sizes 
Lx X Ly. As L —> oo, Ty{L,L + 1) vanishes at a temper- 
ature lower than where Ty{L,L) vanishes. The intersection 
of the dashed line with Ty{L,L) indicates the kink-antikink 
unbinding temperature Tu, — 0.71J predicted by Eq. 1151 . 



we take, 



(16) 



then Ty{kx) is defined bjiSiSi, 



III. NUMERICAL RESULTS 



1 



L^Ly dSAk^dSA^k^ 



(17) 



A. Helicity modulus 

We now present our numerical results. At each tem- 
perature, our simulations consist of typically 10^ — 10^ 
ordinary MC passes through the entire lattice for the 
largest system sizes. In Figs. O and 0] we plot the helic- 
ity moduli Ta; and Ty vs. T, as computed by Eq. © in 
an ensemble with fixed twists A^; = Ay — 0. We show 
an "ordinary" case (no domain wall at T = 0) of size 
64 X 64, which is large enough that any finite size effects 
are negligible for the temperatures shown. In compari- 
son, we also show several "anomalous" cases (percolating 
domain wall at T = 0) of sizes L x {L + 1). For the or- 
dinary case, = Tj, and the bulk transition (where 
the jump discontinuously to zero in the thermody- 
namic limit) is^ at Tkt — 0.81J. In comparison, as L 
increases in the anomalous case, T^; (parallel to the do- 
main wall) in Fig. 121 approaches the value of the ordinary 
case, and so presumably vanishes at the same Tkt- How- 
ever the curves of Ty (transverse to the domain wall) in 
Fig. 2] clearly decrease below that of the ordinary case, 
and presumably vanish in the thermodynamic limit at a 
lower Tuj. 

The reduction seen in Ty for the L x {L + 1) systems, 
as compared to the L x L system, shown in Fig.01 is due 
to the kinks at the Ising-like domain wall. To explicitly 
see this, we can consider the helicity modulus at finite 
wavevector, Ty{kx), defined as the response to a small 
sinusoidal perturbation in the vector potential Ay(vi). If 



In view of the discussion following Eq. 101, equating the 
application of a uniform twist A^ to the addition of a 
constant to the gauge field A^, the helicity modulus Tj, of 
Eq. (|H1) can also be viewed as the zero wavevector helicity 
Ty{kx = 0). In the CG representation, Ty{kx) becomes 
the usual formula for the wavevector dependent inverse 
dielectric functioni^*i^, 



Ty{kx)/J=l 



4^2 J (^q^kx)q{-kx)) 



T 



LxLy fc^ 



(18) 



is the Fourier transform of the 



where g(k) = ^ 
charge distribution. 

Unlike Ty of Eq. ijSJl, which measures the response 
to a uniform twist applied at the boundaries, Ty{kx) 
measures the response to a spatially varying twist ap- 
plied throughout the bulk of the system. For a homo- 
geneous system with periodic boundary conditions, one 
in general expects Ty = limfe^^o Ty(fca;), since the spa- 
tially varying twist becomes uniform as fc^, 0, and 
5A}^_^ — > lS.yl Ly. For free boundary conditions however, 
where the phase angle 0{x, Ly) is not coupled to the phase 
angle 9{x, 0), this equality does not hold. For free bound- 
ary conditions, the absence of any constraint (such as 
in Eq. (jSJ) relating 9{x,Ly) to d{x,0) means that the 
phase angles are free to untwist any additive constant to 
the gauge field, Ay{ri) — > Ay{ri) + Ay/Ly, by choosing 
9{x, y + 1) — 0{x, y) — Ay/Ly; hence if one computes Ty 
by Eq. (jS)! in a free boundary ensemble, one necessar- 
ily has Ty = at any temperature. For the spatially 
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varying twist of Eq. H16() however, no such transforma- 
tion is possible since the perturbing twist is a strictly 
transverse vector function, while the phase angle differ- 
ences give a strictly longitudinal vector function. In this 
case one finds that liuik^^o'^yikx) has the same value, 
as i — > oo, that one has for the system with periodic 
boundary conditions. 

We expect a similar effect to be true in our present 
case. The kink-antikink pairs confined to the one dimen- 
sional Ising-like domain wall can be viewed as a relax- 
ation of the boundary condition. They can unwind, or 
soften the energy of a uniform twist Aj, applied at the 
boundary, but cannot unwind a spatially varying twist 
dAk^ applied throughout the bulk of the system. We 
therefore expect that, as L — » oo, \iiRk^—,QTy(kx) will 
equal the value of Tj, obtained for an ordinary Lx L sys- 
tem, representing the stiffness of the bulk of the system 
on either side of the domain wall; Tj,, however, will be 
a lower value including effects due to the polarization of 
the kink-antikink pairs localized to the domain wall. 

In Fig. we plot, at several different temperatures 
around T^, Ty[kx) for finite kx for L x [L -\- \) systems 
with sizes L = 6A and 128. When we plot the results ver- 
sus the scaled axis of Lkx, we see that the data for the 
two different system sizes collapse to essentially a com- 
mon curve, u{Lkx), at each temperature. As L — > oo, 
such scaling implies that limk^—,o[lmiL^oa^y{kx)] = 
liia^^ao u{k) , is different from Tj^, just as we have ar- 
gued. We may estimate lim^^oo u{k) by taking the max- 
imum value of Ty{kx) for each system size L. In Fig. |^ 
we plot T'^^^ = maxk^[Ty{kx)] versus 1/L for the same 
temperatures as shown in Fig. |S| We give results for 
L = 32,64 and 128. At 1/L = we plot the value of 
Ty{L,L) for L — 64, representing the large L limit for 
the helicity modulus of the ordinary system without the 
Ising domain wall. We see that the values of T™'^'' extrap- 
olate perfectly to Ty{L,L) as 1/L 0. This confirms 
the following conclusion. As L ^ oo in an L x (L -I- 1) 
system, the finite wavevector helicity Ty{kx) is equal to 
the corresponding helicity modulus of an ordinary L x L 
system for all finite values of kx ; this measures the stiff- 
ness of the bulk of the system on either side of the Ising 
domain wall, and is unaffected by the kinks on the wall. 
However the zero wavevector response Tj, to a uniform 
twist Ay , shown in Fig. 01 is softened, and above re- 
duced to zero, by the polarization of the kinks on the 
domain wall. 



B. Finite size dependence 

Returning to Fig. ^ we have indicated the kink- 
antikink unbinding transition temperature T^, as pre- 
dicted in Eq. (|15|l . by the intersection of the line AT/nJ 
with the helicity modulus of the ordinary L x L system. 
This gives an estimate of T^, ~ 0.71J. This value oc- 
curs noticeably above the point where many of the curves 
Ty(L,L + 1) appear to cross. Such a crossing point, if 



1.0 




FIG. 5: Wave vector dependent helicity modulus Tyik^) vs. 
Lkx at several different temperatures T for system sizes L x 
[L -\- 1), with L = 64 (solid symbols), and L = 128 (open 
symbols). The data for the different L collapse to a common 
curve at each T. 




FIG. 6: T;;^''" = maxfcJTy(fc^)] vs. 1/L at several different 
temperatures T. The values at 1/L = correspond to Ty for 
an ordinary 64 x 64 system. Solid lines are fits to a quadratic 
polynomial. 



remaining constant as L increases, is generally taken as 
an estimate for the phase transition in a 2D XY system. 
To explain the difference between this crossing point at 
0.60 J and the above estimate T^, ~ 0.71 J, we need to 
examine the finite size dependence of Tj,(L, L -I- 1) more 
carefully. 

To get the most accurate results, we have found it bet- 
ter to work in the fluctuating twist ensemble and compute 
the phase coherence parameter IS.F of Eq. (O, rather 
than work with periodic boundary conditions and mea- 
sure Ty. In Fig. 13 we plot our results for AL" versus 
system size 1/L, for various temperatures T above and 
below Tw We use L x (L -I- 1) systems with L ranging 
from 32 up to 512. We note that at low T, the behavior 
of A_F is non-monotonic; as L increases, the system first 
softens with a decreasing AF, but then stiffens again as 
AF reaches a minimum and then increases. We denote 
the system size at this minimum by ^fc. Since the sys- 
tem becomes stiffer on length scales L > ^fc, we assume 
that S.k determines the size of the largest kink-antikink 
pairs on the domain wall. At higher T, the minimum 
in AL" disappears, and AF continues to decrease as L 
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increases. The temperature that separates the two be- 
haviors is somewhere between 0.68J and 0.74J, in good 
agreement with the resuh Tu, ~ 0.71J from Fig.Q] based 
on the theoretical prediction of Eq. 115|) . 

As a further check on our resuhs, wc have also directly 
simulated a ID neutral system of logarithmically inter- 
acting charges gkink = ±1/2. The ID interaction poten- 
tial is taken as the 2D L x L lattice Coulomb potential, 
as in Ea. (|ll() . but with height separation fixed a,t y — 0. 
We use an interaction coupling constant of 27rTj,(L,L), 
with Ty(L,L) obtained from our simulations of the or- 
dinary L X L 2D FFXY for L = 64, in order to model 
as closely as possible the interaction between kinks on 
the domain wall in the true L x {L + 1) FFXY system. 
Within this ID simulation we measure the normalized 
histogram of the total dipole moment, P{px), and use it 
to construct^ what would be the free energy F{/S.y) of 
the corresponding L x (L + 1) 2D FFXY system. 



-Tin 



-)/T 



-|-const. 



(19) 

where Vy is the Villain function as in Eq. (|12|l . and 
"const." is a constant term independent of Aj^. 

In Fig. El we plot the resulting AF = F{tt/2) - F{0) 
versus 1/L for the same temperatures and sizes L as in 
Fig- El The agreement between Figs. [7|and|Hlis not ex- 
act, since the coupling between kinks is only equal to 
2TTTy{L, L) on large length scales; the true screening of 
the kink interaction due to charge excitations in the bulk 
on either side of the domain wall is length scale depen- 
dent. Moreover, the domain wall in the 2D FFXY model 
is not a strictly straight one dimensional line; the rough- 
ness of the domain wall (see following section) means that 
height fluctuations can add to the distance of separation 
between kinks. Nevertheless, the agreement is qualita- 
tively very good, indicating again that it is the polariza- 
tion of kink-antikink pairs along the Ising domain wall 
that is responsible for the decrease in the phase stiffness 
transverse to the domain wall. 

One evident feature of both Figs. [7| and El is the very 
large finite size effect. The asymptotic behavior of AF 
only sets in at quite large length scales. Equivalently, the 
correlation length of the kink-antikink pairs, ^k, grows 
large weU below ~ 0.71J. Fitting the data of Fig.Qto 
a quadratic in InL, and determining from the minima 
of these fitted curves, we plot the resulting versus T in 
Fig. El The rapidly growing means that it is difficult to 
get a very precise estimate of T^, directly from the data of 
Fig. \7\ without going to prohibitively large system sizes 
L ^ 512. We believe that this is also the reason that Lee 
alii report a lower value^^ of T/J = 27r(0.09 ± 0.01) = 
0.57 ± 0.06 for the "roughening transition" temperature 
in their dual 2D Coulomb gas. Lee et aZ.'s simulations 
are on a lattice of size L = 64. From Fig. El we see 
that ~ 32 at T/J ~ 0.58, hence kinks already look 
unbound above this temperature for such a small system 
size. The large values of can also be compared to 
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FIG. 7: AF of Eq. (EJ vs. l/L at various temperatures for 
the 2D FFXY model of size L x (L + 1). 
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FIG. 8: AF of Eq. (EJ "^s. l/L at various temperatures for as 
estimated from the dipole histogram of a ID interacting kink 
model (see text). 



other length scales in the system. The correlation length 
of the Ising-like order parameter, gives the typical size 
of an Ising-like domain excitation in an ordinary L x L 
FFXY model. From Ref. (see Fig. 15) we find that ^/ 
is quite small below ~ 0.71J; in particular < 2.5 
for T < 0.77 J. Thus we expect that kink-antikink pairs 
on the domain wall enclosing a typical thermally excited 
Ising-like domain are always effectively unbound. The 
energetics of kink-antikink unbinding will only effect the 
morphology of domains that are much larger than those 
due to typical thermal excitations. 
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FIG. 10: Domain wall width squared, W^, vs. L for several 
different T. Straight lines are linear fits, showing rough do- 
main walls at all T, even below ~ 0.71J. Inset shows the 
slopes of the curves, dW^/dL, vs. T. 



C. Roughening 

Finally, we consider the roughness of the domain wall. 
As noted by KoshunoA*^, an isolated step on a domain 
wall of height two (see Fig. QJi) carries no net charge, 
each corner giving opposite 1/4 charges that therefore 
cancel. Such isolated height of two steps therefore cost 
finite energy, and should roughen the domain wall at any 
finite temperature. To verify this we have explicitly mea- 
sured the domain wall width squared. 



(20) 



where s is an index that counts horizontal length traveled 
along the domain wall (s = a; if there are no overhangs) , 
y{s) is the height of the domain wall at index s, y is 
the average height of the domain wall in the particular 
configuration, and is the number of horizontal links 
in the domain wall (Ns — L ii there are no overhangs). 
To avoid problems with periodic boundary conditions, 
we always measure the domain wall height as relative to 
some initial starting position. In Fig. ^| we plot vs. 



L for various temperatures. The linear growth in as 
L increases indicates that the domain wall is rough for 
all temperatures shown, including T < T^. The domain 
wall diffusion constant, dW'^/dL, shows no observable 
singularity at ~ 0.71J (see inset to Fig. 



IV. DISCUSSION 

Our numerical results demonstrate the existence of the 
unbinding transition for kink-antikink pairs along the do- 
main wall of Ising-like excitations in the 2D FFXY model, 
and that the behavior of this transition is in good agree- 
ment with the theoretical predictions of Korshunov^. We 
show that domain walls are rough at all temperatures, 
and therefore argue that the "roughening transition" 
claimed by Lee and co-worker^iSiii is really the kink- 
antikink unbinding transition. We show that the effects 
of kink-antikink pairs are not readily apparent for Ising- 
like domains such as are typically present due to thermal 
excitation; because of the large length ^fc, such effects are 
important only for much larger domains. 

One of Korshunov's main motivations for investigat- 
ing the kink-antikink unbinding transition was to argue 
that such a transition necessarily implies the existence 
of two separate bulk transitions, Tkt < Tj. We now 
present our own thoughts on this issue. In the original 
paper by Teitel and Jayaprakash'^, two possibilities were 
considered, Tkt < Tj and Tkt ^ Tj. In discussing the 
first case, Teitel and Jayaprakash presented the follow- 
ing scenario. In terms of the dual CG model, the helicity 
modulus T gets reduced from its T = value by fluc- 
tuations that produce dipole moments. If Ising domains 
of typical size ^/ carried a total dipole moment propor- 
tional to their size, they would drive T —^ continuously 
due to the diverging ^/ as T ^ T/ from below. However, 
in addition to these domain excitations, there are also 
pair excitations. The original Kosterlitz-Thouless insta- 
bility criteria would imply that pairs unbind once T falls 
below the critical value T(T) = 2T/tt, which must hap- 
pen at some Tkt below Tj. However, if one computes 
domain energies at T = (and presumably the same 
holds for domain free energies at low T), one finds that 
it is only the domains with vanishing total dipole mo- 
ment that have energies which scale with the perimeter; 
i.e. the only domains which are "Ising-like" are those 
which carry no dipole moment and so cannot give any 
reduction in T! Fortunately, once T increases above Tu,, 
kink-antikink pairs on the boundary of the domain are 
free to unbind, and the Ising domains can now acquire 
large dipole moments at no cost in free energy. The sce- 
nario of Teitel and Jayaprakash is now restored. This 
conclusion is in complete agreement with the numerical 
work of Olsson^, who finds two distinct transitions at 
Tkt < Tj, and argues that the non Ising-like critical be- 
havior claimed by somei at Tj is in fact an artifact of 
finite size effects. 

The kink-antikink unbinding transition may also have 



8 



implications for the non-equilibrium steady state behav- 
ior of the system when the vortices are driven by a 
uniform force, such as is the case for a fully frustrated 
Josephson junction array in an applied uniform d.c. cur- 
rent /. Since / couples linearly to the total dipole mo- 
ment of a domain the force will can can lead to an 
instability^* in domain growth provided the free energy 
of exciting the domain scales less than linearly with p. 
For T < T^^, the binding of kinks to antikinks prevents 
large dipole moments from building up on domains. The 
Ising-like domains, whose free energy scales like the do- 
main length £, are only those domains whose total dipole 
moment vanishes, and hence these domains remain sta- 
ble, and the Ising-like order should persist at small drives. 
For T > Ttu, kink-antikink pairs can unbind, resulting in 
domains whose dipole moment may scale at least propor- 
tional to their length £. In this case, when the current 



exceeds an amount proportional to the Ising domain sur- 
face tension, domains will become unstable to growth and 
the Ising-like order will be destroyed. 
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